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Abstract 

We consider a spherical volume of hot and dense hadronic matter (fireball) expanding into a 
vacuum. It is assumed that initially the fireball matter is in local thermal and chemical equilibrium 
with vanishing collective velocity. The time evolution of the fireball is studied in parallel within 
the GiBUU transport model and an ideal hydrodynamic model. The equation of state of an 
ideal hadronic gas is used in the hydrodynamic calculation. The same set of hadronic species is 
used in transport and fluid-dynamical simulations. Initial coordinates and momenta of hadrons in 
transport simulations have been randomly generated by using the Fermi and Bose distributions 
for (anti)baryons and mesons. The model results for radial profiles of densities and collective 
velocities of different hadronic species are compared at different times. We find that two considered 
models predict essential differences in time evolution of hadron abundances, which are especially 
pronounced for hyperonic species. This gives an evidence of a strong deviation from chemical 
equilibrium in expanding hadronic matter. 

PACS numbers: 24.10.Lx, 24.10.Nz, 25.75.Dw, 25.75.Ld 
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I. INTRODUCTION 



Physics of high energy heavy-ion collisions exhibits a rapid development during last 
several decades. A lot of efforts has been made to extract the information on properties of 
hot and dense nuclear matter from experimental data. Numerous theoretical models have 
been proposed to describe the complicated dynamics of multi-particle systems produced in 
nuclear collisions. An important place in understanding the main features of experimental 
data is still occupied by simple fireball jl, 2], blast wave 3, 4| and thermal 0, 6| models. 

More sophisticated hydrodynamic models (see e.g. a recent review have been very 
successful in reproducing a large amount of experimental data in a wide range of bombarding 
energies. An attractive feature of these models is their capability to study the sensitivity of 
observables to the equation of state (EoS) of strongly interacting mater and, in particular, 
to a possible deconfinement phase transition. The hydrodynamic models (HDM) are usually 
applied to describe the evolution of matter only at some intermediate stage of a heavy-ion 
collision. They are not well suited for early and late stages since large deviations from local 
thermodynamic equilibrium are expected in this case. 

As a rule, these models assume the formation of a quasi-equilibrated fireball at an inter- 
mediate stage of the collision process. The geometrical and thermodynamic parameters of 
the initial fireball are normally chosen to achieve the best fit of observed data. The (2+1)- 
dimensional HDM js-12] became especially popular in recent years. These models describe 
the dynamics of a cylindrical fireball expanding in transverse and longitudinal directions. 

Of course, the hydrodynamic approach can not be directly applied to simulate the late 
stages of a heavy-ion reaction when collisions of particles become too rare to maintain 
the thermodynamic equilibrium. The assumption of an instantaneous transition to the 
collisionless propagation of particles (" freeze-out" ) has been introduced [3, [l^| to calculate 
the asymptotic particle spectra in HDM. However, the direct kinetic simulations of space- 
time distributions of last-collision points in relativistic nuclear reactions show [lij] that the 
freeze-out is in fact a continuous process, even in central interactions of heaviest nuclei (for 



details, see 0, Q). It has been also conjectured [17|, |18| that the stages of a "chemical" 
(saturation of particle multiplicities) and "kinetic" (saturation of particle spectra) freeze- 
out should be separated in time. It is believed that the chemical freeze-out takes place 
earlier because of much faster drop of inelastic collision rates as compared to elastic ones 
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in expanding matter. On the other hand, there are several fluid-dynamical studies [19|, |20 | 
where experimental particle spectra have been successfully reproduced without introducing 
such a separation. 

Thermodynamic equilibrium is not postulated in microscopic transport models of nu- 



clear collisions. Currently, the UrQMD 21], GiBUU [22J, |23|, HSD [2J, |25|, QGSM [26] 



and AMPT 



271 ] transport models are widely used for simulating relativistic nuclear col- 



lisions. These models are not well suited for describing high density states of hadronic 
matter, where multiparticle (non-binary) channels of hadronic interactions are presumably 
important. As a rule, contributions of such channels are disregarded in existing versions of 
transport models. On the other hand, as demonstrated in Refs. 28N30|. multimesonic chan- 



nels of baryon-antibaryon pair production might be responsible for high yields of antiprotons 
and antihyperons observed at SPS and RHIC energies. 

A more realistic description of the freeze-out pro cesses in heavy-ion collisions can be 



achieved in a hybrid "hydro-cascade" model 11 
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33| . In this approach, hydrody- 



namic and cascade simulations are applied, respectively, for intermediate and late stages 
of the reaction. The hybrid model implicitly assumes the existence of a space-time region 
where the hydrodynamic and cascade simulations give similar results. The characteristics 
of hydrodynamic flow taken at a certain hypersurface are used for generating phase-space 
coordinates of hadrons which serve as an input for subsequent transport simulations. 

Although such a procedure seems to be intuitively justified, it should be verified for more 
or less typical situations. In the present work, a spherical expansion of an ideal hadronic gas 
is simulated in parallel by using the transport and hydrodynamic models. Specifically we use 
a non- viscous hydrodynamics with an ideal-gas EoS and the transport GiBUU model 
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23] 



(without mean-field interactions). The initial fireball is assumed to be in thermodynamic 
equilibrium. Thermodynamic parameters of this state are used to generate randomly the 
set of initial hadrons for subsequent transport calculations. By the direct comparison of 
results predicted by the GiBUU and HDM we investigate possible deviations from thermal 
and chemical equilibrium in the course of expansion. 



34-; 



which studied the particle 



Finally, we would like to mention several earlier works 
composition in an expanding fireball within a " hadrochemical" approach. There a set of rate 
equations for particle abundances was combined with a simplified hydrodynamic description. 
But the calculations had been performed for a spatially uniform background neglecting 
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surface effects. The direct comparison between hydrodynamic and transport calculations 



37] 



have been recently made in the case of a cylindrical gluonic fireball in Ref. 

The paper is organized as follows: parameters of the initial fireball are given in Sec. IIA, 
our hydrodynamic model is formulated in Sec. IIB. In this section we also show radial 
profiles of densities and collective velocities obtained by numerical solution of fluid-dynamical 
equations. Our transport model is formulated in Sec. IIC In Sec. Ill we compare the 
predictions of the hydrodynamic and transport models. Our conclusions are presented in 
Sec. IV. The procedure of initial event generation is described in Appendix. 

II. DESCRIPTION OF THE MODELS 
A. Initial state 



We assume that at t = the fireball is locally equilibrated and has vanishing collective 
velocity. The initial radial profiles of energy (e) and baryon (n) densities are parametrized 
by the Woods-Saxon distribution with some radius R and diffuseness a: 

(1) 



r-R , 
exp | ) + 1 



a 



e (r,0) = e W(r), n(r, 0) = n W(r), W(r) 
Below we choose the following parameters: 

R = 6 fin, a = 0.3 fm, e„ = 1.3 GeV/fm 3 , n = 0.45 fm" 



(2) 



The initial values of energy and baryon densities at the fireball center, and no, are typical 
for high density states of hadronic matter created in heavy-ion collisions at — 10 AGeV 
(see e.g. Ref. |38j). 



B. Hydrodynamic simulation of fireball expansion 



In this section we describe the HDM for simulating the fireball expansion. The model as- 
sumes that deviations from local equilibrium are small at any space-time point (r, t). In this 
approximation, the single-particle phase-space distribution function (DF) of 2-th hadronic 
species is equal to a locally equilibrated DF characterized by certain temperature T, 
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chemical potential fa and collective 3- velocity v. Considering the hadronic system as a mix- 
ture of ideal gases of (anti)baryons and mesons, one can approximate f- by the Fermi or 
Bose DF (h = c= 1): 



(2tt)* 



Ei 



exp 



fa 



T 



± 1 



(3) 



where gi is the spin-isospin degeneracy of z-th hadronic species |54j, Ei = j(Ei — pv) is 
a single-particle energy in the local rest frame (LRF), Ei = \frn% + p 2 is the corresponding 
energy in an arbitrary frame, m, is the mass of i-th hadrons, and 7 = (1 — v 2 )^ 1 ^ 2 . Plus or 
minus in the r.h.s. of Eq. ([3]) correspond, respectively, to fermions or bosons. 

In our fluid-dynamical calculations we assume a fireball matter to be in chemical equi- 
librium with respect to strong interactions and decays of hadrons. In this case, fa may be 
represented by linear combinations of the baryon (fa) and strange (fas) chemical potentials: 



fa = BijjL + Si fas 



(4) 



where Bi = 0, ±1 and Si = 0, ±1, ±2 . . . are, respectively, the baryon and strangeness quan- 
tum numbers of the i-th hadrons. 

The values of the net baryon (n), strangeness (rig) and energy (e) densities in LRF, as 
well as pressure (P) may be expressed [39| in terms of integrals of the DF ([3]) over the 
3-momentum in this frame: 



(n \ 

e 

\P ) 



Si 

Ei 
P 2 



(5) 



\ 3Ei/ 

where the sum is taken over all hadronic species. Assuming that the strangeness den- 
sity n s = at any space-time point, and using Eqs. (J3]) (J5j) , one can express all thermo- 
dynamic quantities (e.g. T, fa fas, P) in terms of two independent variables n and e. In 
particular, in this way we obtain the EoS of the fireball matter, P = P(n,s), which is used 
in our HDM (see below). 

To facilitate the comparison with transport simulations, in our hydrodynamic calculations 
we use the EoS with the same set of hadrons and hadronic resonances as in the GiBUU 

0) 



model (see Appendix A of Ref. 23]). The set of nonstrange baryons (Bi 



I, Si 



TABLE I: Hadronic species included in the calculation. 
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includes nucleons (iV), isobars (A(1232)) and their excited states. We also take into account 
hyperons stable with respect to strong decays (Y = A, S, H, Q~) and hyperonic resonances 
(Y* = S(1385), A(1405), . . .). In addition to baryons we include corresponding antibaryons 
(B{ = —1). The mesonic set (Bi = 0) consists of stable mesons ir, K, K, 7], r\' and resonances 
p,u,a,(p,K*,K*. 

The list of baryons and mesons included in our calculations is given in Table d 55 ]. 
The thermodynamic properties of a fireball matter, in particular, its EoS, are obtained 
assuming zero widths of resonances. Using Eqs. (EE))-© and characteristics of hadrons 
from Table H we calculate the initial radial profiles of temperature T, chemical poten- 
tials /j., Us and partial densities of hadronic species, rij = J d 3 pf} . This calculation 
gives the values T ~ 180 MeV,/i ~ 328 MeV, fj, s ^ HO MeV at the fireball center (r = 0). 
The last two columns of Table [J show the equilibrium multiplicities of i-th hadrons (iVj) 
and antibaryons (Nj) in the initial state. These quantities are obtained by the vol- 
ume integration of the corresponding partial densities. One can see that most abundant 
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hadrons at t = are pions, nucleons and isobars. The total multiplicity of initial an- 
tibaryons (approximately 18.6) is much smaller than the total baryon number of the fireball 
#tot = E \Bi\(Ni -Nj) = j c/ 3 rn 7 ~ 417. 

i 

To describe the dynamics of spherical fireball expansion, we solve the equations of the 



^l + l)-dimensional ideal hydrodynamics. They can be written in the form 40]: 



dt \dr r 

£ + (» + H) ( .* + „-o. 

Here v = vr/r is the radial component of fluid velocity, M = 717 is the net baryon density 
in the fireball cm. frame, E and M are the components of the energy-momentum tensor: 

E = T 00 = 1 2 (e + v 2 P), M = T 0r — v(E + P) . (9) 

The initial conditions for these equations are given by the relations (pQ) and v(r, 0) = 0. 
Solving ([6])-((9]) together with the EoS P = P(n,e) gives the radial profiles of n,e,v at 



fixed t > 0. Similar hydrodynamic studies o 
EoSs have been performed earlier in Refs. 

The numerical solution of Eqs. 
transport algorithm SHASTA 



a spherical fireball expansion by using simplified 



4lJ-|43| 
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IS]) - © has been obtained by using the flux- corrected 
451 ] . Typically, we choose the cell sizes Ar = 2.5 • 10 -2 fm, 



At = 5 ■ 10~ 3 fm/c of the (r,t) grid. Similarly to Ref. 38J, in this calculation we use linear 
interpolations of the P(n,e) table prepared with fixed steps in n and e. We have checked 
that our numerical scheme conserves the total baryon number, energy and entropy of the 
fireball with relative accuracy better that 1%. 

The radial profiles of fluid- dynamical quantities calculated for different times t within the 
HDM are shown in Figs. HH21 In particular, thick and thin lines in the left panel of Fig. [1] 
show the profiles of n and A/", respectively. As one can see in Fig. [TJ a rarefaction wave 
propagates from the fireball periphery and reaches the fireball center at t ~ 10 fm/c. At 
later times the matter in the whole fireball is involved in the expansion. At such times the 
radial profile of collective velocity can be approximated as v = Hr a with a < 1. In Fig. [2], 
we do not show external parts of velocity profiles corresponding to dilute regions of matter 
with n < 10~ 3 fm -3 . These parts certainly can not be realistically described within the ideal 
HDM because of large local Knudsen numbers. 
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FIG. 2: The radial profiles of the collective velocity (left panel) and temperature (right panel) 
calculated in HDM. The outer parts of velocity profiles corresponding to densities n < 10~ 3 fm~ 3 
are omitted. The dots in the left panel correspond to spatial points where n = 10~ 2 fm~ 3 . 
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C. The GiBUU transport model 



To study possible deviations from local equilibrium, the expansion of the hadronic fire- 
ball has been also simulated within the GiBUU transport model (version 1.2.2) 22l |23|. |46|. 



The detailed description of the model is given in Ref. [23|, below we present only its short 
summary. We apply GiBUU in a cascade mode, i.e. omitting contributions of nuclear and 
Coulomb potentials. The evolution of hadronic phase-space DFs fi(r, p, t) due to free prop- 
agation of hadrons, as well as due to their two-body collisions and resonance decays is de- 
scribed by kinetic equations with corresponding collision terms. They are solved by using the 
standard test-particle method which includes the Hamiltonian equations of motions for 
test particles and their two-body collisions (or decays) generated by a Monte Carlo method. 
In present calculations we use one test particle per hadron and consider a large number of 
parallel ensembles which we call "events". Each event is characterized by a specific sets of 
hadronic species, as well as their coordinates Vj(t) and momenta Pj(t). 
The DF of i-th hadronic species is represented as follows 

fi(r,P, t) = (j^S[r- rj (t)\ 5 [p - Pj (t)]J , (10) 
where the angular brackets denote averaging over all events and the sum runs over all 



hadrons of the type i existing at given time t in a given event 56j. The hadrons propagate 



along straight-line trajectories (rj — pj/Ej, pj = 0) between their two-body scatterings. 



The latter change the momenta abruptly and also lead to production of new hadrons 57]. 
Resonances are allowed to decay during time evolution. 

For calculating collective velocities of particles, one should know radial profiles of the 
4-current density Jf and the energy- momentum tensor T? v of the hadronic species i. As 
in the HDM, we consider initial states with spherically symmetric distributions of particles. 
Let us consider a thin shell occupying a region of points with radii between r and r + Ar. In 
the spherical coordinate system with the origin at the fireball center, only the components of 
Jf , T? v with /i, v = 0, r are nonzero in the limit of a large number of events. By using Eq. ( ITOl) 
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one can get the relations: 

r I N * \ 

(11) 
(12) 

kj=i / 

where Qj = (47rr 2 Ar) _1 if the radial coordinate r 3 - of hadron j in a given event falls into the 
interval (r, r + Ar) and Oj = in the opposite case. The second equalities in (ITT]) — (fl2il are 
obtained in the limit of small Ar and large number of events. A similar expression for J[ 
takes place after replacing Qj — > fjO^ in (ITT]) , where is the radial component of the j-th 
hadron velocity. The replacements of Qj by vJQj and (vj) 2 0j in (I12p give, respectively, the 
corresponding relations for T® r and T[ r . 



For generating initial events we use a procedure similar to that suggested in Ref . [33( . It is 
described in Appendix. These events have been generated according to the Fermi and Bose 
DF of hadrons ([3]). In principle, in this way one should obtain the same ensemble averaged 
initial distributions in both models (58 



Unless otherwise stated, the GiBUU simulations 
are averaged over 5000 events. Note that in our generating procedure we do not fix the 
total strangeness and charge of particles in a single event. However, we have checked that 
fluctuations of these quantities are relatively small for typical events. In particular, absolute 
values of net total strangeness do not exceed 1% of B tot . 



III. RESULTS 

A. Comparison of density profiles 

The process of fireball expansion into a vacuum has been simulated in parallel within 
the hydrodynamic and GiBUU models. In this section we compare the density profiles for 
different hadronic species. Figure |3] shows the nucleon and pion densities in the fireball 
cm. frame. These densities are summed over different isospin states of corresponding 
hadrons. Histograms in Fig. [3] are obtained within the GiBUU model by using Eq. ffTT]) 
for i = N, 7r. Thin lines represent the profiles of densities obtained within the HDM. 
In the case of pions, deviations between two calculations are visible already at relatively 
short times t < 5 fm/c. In particular, one can clearly see an excess of pions predicted 
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by GiBUU in a central region. This effect can be explained by multi-hadron absorption 
processes missed in GiBUU but implicitly included in the HDM. Another possible reason is 
the Bose-enhancement effects neglected in transport calculations, but taken into account in 
generating initial events (see the footnote at the end of Sec. Ill Cp . At later times the HDM 
densities are systematically larger than the GiBUU predictions in this region. 




5 10 15 20 5 10 15 20 

r (fm) r (fm) 

FIG. 3: The radial profiles of nucleon (left panel) and pion (right panel) densities calculated 
within the GiBUU (histograms) and hydrodynamic (thin lines) models. 

The density profiles of A hyperons and kaons are shown in Fig. HI One can see again that 
in contrast to the HDM, the particle densities in the central region increase in GiBUU at 
early times. However, at t > 10 fm/c, the local densities of A's predicted by HDM decrease 
much faster with time than those obtained in GiBUU. This is a consequence of the rapid 
drop of the activation factor exp [(/ia — m\)/T] as the matter cools down. This factor is 
especially important for massive particles like hyperons and baryon resonances. As we shall 
see in Sect. IIII Cl a faster drop of densities gives rise to a qualitatively different time evolution 
of heavy particles' multiplicities in the HDM and GiBUU. 

B. Collective velocities 

There is no unique way to determine a collective velocity (CV) of particles in transport 
models. One possibility is to define this velocity (we denote it by vi) as the velocity of LRF 
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FIG. 4: Same as Fig. [3l but for radial profiles of A hyperons (left panel) and kaons (right panel). 

where the 3-vector of the current density vanishes in average over many events. In the case of 
spherically symmetric expansion, the CV has only the radial component. In general, within 
GiBUU, the hadrons of different types have different CVs. Using the Lorentz transformation 
to the frame where the radial current of i-th hadronic species equals to zero, one obtains 
the relations: 



(i) 



J? 



(13) 



The current density components in second equation are calculated by using ( TTTj) . 

Another possibility is to find the velocity of reference frame with vanishing 3-momentum 
flux. We denote the corresponding CV by Vi- From the Lorentz transformation for the 
energy-momentum tensor one has 

2 ^ Or 



i0r 



-> v 



(14) 



It can be shown that in the case of a local equilibrium both definitions give the same values 



of CV: v\ 



(i) 



v, where v is the flow velocity in the HDM. Our GiBUU calculations 



show that Vi and ?4 are nearly equal for all kinds of hadrons. In the following we determine 
the CV by using the definition (TTBl) . 

The radial profiles of CV for nucleons and pions are shown in Fig. [5] for several times t. 
One can see that these CV are not equal to each other in the transport calculation. At 
small t significant differences between and the hydrodynamic velocity exist only at 
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FIG. 5: Collective velocities of nucleons (left panel) and pions (right panel) as functions of radius r 
at different times indicated in the left panel. Thick and thin lines are calculated within the GiBUU 
and HDM, respectively. Outer parts of profiles corresponding to densities below 10 -3 fm -3 are 
omitted. 

the fireball periphery. At later times the differences become visible for all r. The deviations 
between the models are especially large for nucleons at t > 10 fm/c. These deviations can be 
explained by dissipative effects effectively included in GiBUU, but disregarded in the HDM. 



C. Hadron abundances 



In this section we compare evolutions of particle multiplicities predicted by the GiBUU 
and HDM. These multiplicities are obtained from radial profiles of partial densities by inte- 
gration over the fireball volume. In addition to the partial densities of "free" hadrons, rij, 
below we also calculate the corresponding total densities, n\ ot , which include hadrons "hid- 
den" in heavier hadronic resonances 



13]: 



tot 



11, 



(15) 



where is the average number of i-th hadrons, produced in decays j — > i + X. The sum 
in the r.h.s. runs over all resonances included in our simulations (see Table [I]) and having 
strong decays into i-th hadrons. The coefficients d* are calculated using branching ratios 
given in Ref. 481 ] . 
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FIG. 6: Multiplicities of pions, kaons and antikaons as functions of time in the hadronic fireball 
expanding into vacuum. Thick (thin) lines are calculated within the GiBUU (hydrodynamic) 
model. 7 r tot ) jf tot and R tot are total numbers of pions, kaons and antikaons including mesons, 
hidden in hadronic resonances (see Eq. (|15p ). 
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FIG. 7: Same as Fig. El but for nucleons, A-isobars and antibaryons. _/V tot is the total number of 
nucleons including those, hidden in hadronic resonances. 
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FIG. 8: Same as Fig. [7J but for multiplicities of hyperons. Y(Y*) denote multiplicities of all 
stable (unstable) hyperons. A tot is the total number of A's including those, hidden in hyperonic 
resonances. 
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FIG. 9: Total production (dashed line) and absorption (dash-dotted line) rates of A hyperons as 
functions of time in GiBUU. The solid line shows the net (production minus absorbtion) rate. 



Figures [6H8] show the time evolution of hadron multiplicities predicted by the GiBUU 
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and hydro dynamic models. One can see a qualitative difference between the results of two 
calculations. In particular, the GiBUU model predicts nearly constant total multiplicities of 
pions, kaons and hyperons including hadrons hidden in resonances. However, these quantities 
decrease noticeably within the HDM. 

The two models predict very different time evolution of multiplicities of hyperons and 
hyperonic resonances. As one can see from Fig. [BJ the multiplicity of A hyperons noticeably 
increases in GiBUU, but decreases in the HDM. According to our transport calculation, 
the number of A's increases from 13 at t = to about 50 at t = 40 fm/c. Within this time 
interval many A particles are still "bound" in hyperonic resonances. When time increases 
from to 40 fm/c, the relative fraction of bound A's decreases from about 80% to 30% 59]. 
A long duration of A production in GiBUU is directly connected with a slow decrease of the 
abundances of hyperonic resonances (compare dash-dotted curves in Fig. [8j. This in turn 
follows from significant regeneration of these resonances in baryon-baryon and baryon-meson 
collisions. 

0.4 



03 
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20 



FIG. 10: (Anti)kaon to pion ratios as functions of time in expanding fireball. Thick (thin) lines 
are calculated within the GiBUU (hydrodynamic) model. 

To analyze relative importance of different channels of A production and absorption we 
have calculated total rates of such processes in GiBUU. This analysis shows that for our 
parameters of initial state, the most important production channels are decays Y* — > AM 



1(3 



(here M denotes a nonstrange meson) and £5 — > AX, Y*B — > AX reactions (B is a non- 
strange baryon). On the other hand, most significant channels of A absorption are the 
reactions AM — >■ Y*X and AB — > (E,Y*)X. As one can see in Fig. EH both the A produc- 
tion and absorption rates are rather high at t < 20 fm/c, but they nearly compensate each 
other: the net production rate is of the order of 1 c/fm at t < 40 fm/c. Only later the mul- 
tiplicity of A hyperons saturates. It is interesting that large characteristic times (exceeding 
approximately 40 fm/c) of A production have been obtained earlier in the UrQMD llj and 



BUU 



491 ] calculations. 



Figure [TU] shows the time dependence of (anti)kaon to pion multiplicity ratios. Again, 
one can see significant differences between the model predictions. For example, the ratios 
K tot /TT tot and K /vr tot are practically constant in the GiBUU calculation. However, they 
rapidly decrease with time in the HDM. These results may imply large deviations from 
chemical equilibrium in expanding baryon-rich hadronic matter created at FAIR, NICA and 
lowest RHIC energies. 



D. Evolution of energy distributions 




0.2 0.4 0.6 0.8 1.0 1.2 0.2 0.4 0.6 0.8 1.0 1.2 
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FIG. 11: Kinetic energy distributions of pions (left panel) and nucleons (right panel) in expanding 
fireball. Lines show the GiBUU results at different time moments t. 

It is instructive to see, how particle momentum spectra evolve with time in GiBUU. 
In Figs. riTHT2l . we show the kinetic energy distributions of pions, nucleons, kaons 
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and A hyperons at different t. One can clearly see the importance of final state rescatter- 
ings and resonance decays in formation of these spectra. The collective expansion results in 
concave and convex shapes of momentum distributions for pions and nucleons, respectively. 
These effects are less pronounced for kaons and A's. The GiBUU simulations predict the 
sequential behavior of freeze-out. Indeed, the formation of asymptotic distributions is taking 
place at different times for different hadrons. First, this happens for kaons (at t ~ 20 fm/c), 




E kin (GeV) E kin (GeV) 

FIG. 12: Same as Fig. 1111 but for kaons (left panel) and A's (right panel). 



then for pions and nucleons (t > 30 fm/c) and finally, for A hyperons (t > 50 fm/c). One can 
not clearly distinguish the stages of "chemical" and "thermal" freeze-out, especially for A's. 

IV. CONCLUSIONS 

In this paper we have compared the results of transport (GiBUU) and hydrodynamic 
calculations for the expansion of a baryon-rich hadronic fireball. The initial events for 
transport calculations have been generated in accordance with locally equilibrated phase- 
space distributions of hadrons. We have found significant differences in space-time evolution 
of partial densities and collective velocities of hadrons predicted by these models. Also, the 
two models predict very different time evolution of hadron abundances. 

The present study demonstrates that a simple picture of sequential chemical and ki- 
netic freeze-outs is strongly distorted by the decay and regeneration of hadronic resonances. 
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According to our transport calculations, an early saturation of the particle multiplicities 
in an expanding fireball occurs only for the total yields, which include hadrons hidden in 
heavier resonances. This is especially important for strange hadrons. The hydrodynamic 
simulations with chemically equilibrated EoSs strongly underestimate total multiplicities of 
hyperons and kaons. We believe that strong deviations from chemical equilibrium predicted 
by transport calculations for such hadrons are real and they should be explicitly taken into 



account in future hydro-cascade calculations 



60]. 



However, one should be very careful in selecting and combing the corresponding models. 
Obviously, a smooth matching of the models is only possible if there exists a stage when both 
transport and hydrodynamic models predict a similar behavior of macroscopic quantities. On 
the basis of our present analysis we conclude that these models predict an essentially different 
evolution of particle densities and velocities. Therefore, the combined model predictions will 
strongly depend on the choice of the transition hypersurface. 

One should bear in mind, that in this paper we have disregarded viscosity effects in 
the hydrodynamic calculations. These effects may be responsible for the differences in col- 
lective velocities predicted by the two models. In the future we are going to check sensitivity 
of the results to the choice of initial conditions. In particular, we plan to extend our analysis 
to the case of nonzero initial collective velocity of the fireball. It would be also interesting to 
perform calculations with inclusion of relativistic mean fields. In this way one can effectively 
account for soft non-binary hadronic interactions disregarded in our present calculations. 
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Appendix A: Generation of initial events 

In our calculations we assume that initial particles are located inside a cube with dimen- 
sions | a; |, | ?/ 1, |^| < L (r = corresponds to the fireball center). This cube is divided into 
smaller cells with same sizes A/ along each coordinate axis. Below we choose L = 10 fm, 
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Al = 0.2 fm. For each cell we determine average numbers of the i-ih hadrons, Ni = rii (A/) 3 , 
where rij is their partial density at the cell center. The latter is calculated by integrating 
the equilibrium DF fl3]) over the 3-momentum. This calculation is performed for all mesons, 
baryons and antibaryons listed in Table [B 

It is assumed that for each hadronic type particle multiplicities in a given cell are dis- 
tributed in accordance with the Poisson distribution. In this case the probability to find n 
hadrons of the type % in the cell is given by the expression 

-jryn 

«,« = _L<«p (A.l) 

From Eq. ( 1A.1I) one can see, that the probability that a given cell is empty, i.e. contains 
no hadrons, equals P = Y\wq = ex P( — -^ tot ), where N tot = ^2 Ni is the total average 

i i 

multiplicity of hadrons of all types in the cell. For our choice of parameters, ~ < 1 
even for most abundant hadrons. Probabilities to find more than one hadron (n ^ 2) are 
negligible for all cells. In our calculations we consider only events with cells containing no 
more than one hadron. 

At the first step of the initialization procedure we generate randomly particle's coordinates 
x, y, z assuming that each coordinate is homogeneously distributed in the interval (-L, L). 
Then we find the cell containing the point (x,y,z) and decide whether it is empty or not. 
To determine this, we choose a random number £ with the homogeneous distribution in the 
interval [0, 1]. If the inequality £ < 1 — Pq does not hold, the cell is considered empty (it 
is excluded from further treatment of a given event). In the opposite case we consider the 



cell as "filled" and generate the type of hadron contained in it [611 ] . Following Ref. 33] 
we assume that the relative probability of z-th species equals Ni/N tot . For hadrons with 
nonzero isospin / we randomly generate its isospin projection I 3 . This is done assuming the 
homogeneous distribution of /3 in the interval 1/3] ^ /. 

At the second step we generate a hadronic 3-momentum p. It is postulated that the mo- 
mentum distribution of the i-th hadron in a given cell is proportional to the equilibrium DF 
introduced in Eq. ([3]). Normalizing this distribution to unity, one gets the relation 

d 3 w- f {eq) 

^P- = ^— = i Pi (p). (A.2) 
Grp rii 

To generate particle's momentum in accordance with distribution (1A.2I) . we use the well- 
known rejection method 53| extending it to the case of three dimensions. First, we randomly 
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choose three components of p, assuming that each component is homogeneously distributed 
in the interval [— Pm,Pm] where p m = 3 GeV/c. Generating a new random number £ e [0, 1] , 
we accept the momentum p if the inequality £ < <Pi(p) / <Pi(0) holds [62J. 

This procedure is continued for other cells until the total net baryon number of generated 
particles does not exceed the value B tot determined by the initial density profile n(r, 0) 
(see Eq. ([]])). The resulting set of coordinates, momenta and types of particles is considered 
as a single event. It is used as an initial condition for subsequent GiBUU simulations. 
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